Unconventional metallic conduction in two-dimensional Hubbard- Wigner lattices 
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The interplay between long-range and local Coulomb repulsion in strongly interacting electron 
systems is explored through a two-dimensional Hubbard- Wigner model. An unconventional metallic 
state is found in which collective low-energy excitations characteristic of the Wigner crystal induce 
a flow of electrical current despite the absence of one-electron spectral weight at the Fermi surface. 
Photoemission experiments on certain quarter-filled layered molecular crystals should observe a gap 
in the excitation spectrum whereas optical spectroscopy should find a finite Drude weight indicating 
metallic behavior. 

PACS numbers: 71.30.+h; 71.27. +a; 71.10.Fd 
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I. INTRODUCTION 

Coulomb interactions in two dimensions (2D) can 
lead to different forms of electron localization rang- 
ing from Mott insulators 1 - to Wigner crystals 2 -. Half- 
filled narrow-band systems such as cuprate 3 ^ and or- 
ganic superconductors 5 ^ display insulating states in their 
phase diagrams which are not described by band theory 
approaches. This is due to the strong local Coulomb 
repulsion driving the system to a Mott metal-insulator 
transition (MIT) which localizes electrons at each site 
of the lattice. On the other hand, a gas of electrons 
in a positive uniform background can localize through 
Wigner crystallization when the long-range Coulomb re- 
pulsion overcomes the kinetic energy, which is relevant 
to the MIT observed in two-dimensional electron gases 
(2DEG) at sufficiently low densities^. Besides, there are 
many 2D crystals with non-integer filled narrow bands 
which arc neither Mott nor Wigner insulators but rather 
display "Wigner crystallization on an underlying lat- 
tice" (WL)£, with both aspects of electron localization. 
The family of quarter-filled layered organic materials, 0- 
(BEDT-TTF) 2 MM'(SCN) 4 (M=Rb, Cs, Tl,M'=Zn, Co) 
represent clean realizations of WL displaying a subtle 
competition between charge ordering, superconductiv- 
ity and unconventional metallic phases in their phase 
diagram-^. These systems pose the challenging question 
of whether the interplay between short (Mott) and long- 
range Coulomb interactions (Wigner) 11 : 12 on a lattice can 
lead to new ground states and excitations and whether 
such novel states of matter can be experimentally tested. 

Extensions of the Hubbard model including the long- 
range Coulomb interaction have mostly been limited to 
one-dimensional Wigner lattices ^ 13 : 14 ' 15 : 16 . In these sys- 
tems the degeneracy of the classical ground state at non- 
integer fillings leads, in the strongly interacting regime, 
to low-energy excitations consisting of domain walls with 
fractional charge. This type of collective excitations, 



which play a major role in the melting of the Wigner lat- 
tice on the one-dimensional chai n 14 : 15 , have been much 
less explored in two dimensions! 17 : 18 Here we find that 
collective excitations in 2D can give rise to metallic con- 
duction in a charge ordered metallic state (COM) even 
though there is a finite charge gap in the one-electron 
spectrum. By combining conductivity and photoemis- 
sion/tunncling experiments, such an anomalous metallic 
phase could be clearly identified and discerned from a 
more conventional charge density wave (CDW) metal, in 
which the conducting behaviour is ensured by portions 
of the Fermi surface that remain gaplcss. 

The paper is organized as follows. The model and the 
numerical method that we use are presented in Sec. II. 
The main results are contained in Sec. Ill, and their 
possible relevance to experimental systems is discussed 
in Sec. IV. In the appendices we introduce an analytical 
model valid in the limit of strong long-range Coulomb 
interactions, that provides a clear interpretation of the 
numerical results in terms of defects of the checkerboard 
ordering. 



II. MODEL AND METHOD 

The minimal model to describe the effect of the long- 
range part of the Coulomb interaction in 2D is the 
Hubbard- Wigner model (HWM) at one-quarter filling on 
a square lattice: 
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where c ia creates an electron with spin a on site 



(1) 
the 



kinetic energy is parametrized by the hopping amplitude 
t between neighboring sites and the onsitc Coulomb in- 
teraction by U. The long-range contribution is: Vij = 
V/\i — j\, with \i — j\ the distance between two different 
sites on the lattice and V a parameter controlling the 
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FIG. 1: Sketched phase diagram of the quarter-filled 
Hubbard- Wigner model. The uniform metallic phase is sta- 
bilized by the inclusion of long-range Coulomb interactions 
(HWM) compared to the model with short range interactions 
(EHM) and an intermediate charge ordered metallic phase 
emerges in between the uniform metal and the Wigner lat- 
tice. 



strength of the non-local Coulomb repulsion. The model 
is solved using Lanczos diagonalization on finite clusters 
of up to L — 18 sites with periodic boundary conditions. 
The Coulomb potential between two different sites in- 
side the cluster is calculated from Ewald summations 
which account for an infinite periodic array of simula- 
tion cells. In order to single out the effects introduced 
by the long-range part of the interaction we compare re- 
sults of the full model with the extended Hubbard model 
(EHM) containing U and a nearest-neighbor V only. We 
deliberately choose a large U = 100£ in order to disen- 
tangle the relevant low-energy excitations of the system, 
focusing primarily on the charge sector. Nevertheless, 
the main conclusions obtained here remain valid as long 
as U is larger than the bandwidth, 8t, which includes 
the realistic case U « 20t relevant to the organic cystals: 
6>-(BEDT-TTF) 2 X. 
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FIG. 2: Finite-size scaling of the static charge correlation 
function evaluated at (7T, 7r) (a) and Drude weight (b) for 
different values of V. The charge ordering transition is es- 
timated to occur at Vco ^ 4i from the linear extrapolation 
to the thermodynamic limit of C(7r,7r) vs. l/L. The Drude 
weight is plotted as a function of 1/ y/Z, displaying insulating 
behavior for V > Vmi ~ 8t. 



III. RESULTS 

A. Phase diagram 

The phase diagram of the HWM obtained from exact 
diagonalization (ED) of finite-size clusters is illustrated 
in Fig. [U where two distinct phase transitions can be 
identified. The first one corresponds to the charge order- 
ing instability of the homogeneous metal, indicated by 
the presence of non-zero charge correlations at the crit- 
ical wavevector Q = (ir,ir). In Fig. [5] (a) we show the 
finite size scaling of the charge correlation function on 
L = 8, 10, 16, 18 clusters as a function of l/L. The ex- 
trapolation to the thermodynamic limit suggests a charge 
ordering transition at about Vco ~ (4 — 5)t in good 
agreement with previous results using the Path Inte- 
gral Renormalization Group (PIRG) on 4 x 4, 8 x 8 and 
12 x 12 clusters (see Fig. [5] and Ref. HJ. This value 
is much larger than the 1.5t obtained for the EHM (see 
also Ref. Il9f ). This indicates that the uniform metallic 
phase is notably stabilized by the inclusion of long-range 
interaction al 20 . To trace back the physical origin of 
this phenomenon, we observe that the CO instability can 



be understood through the divergence of the charge sus- 
ceptibility in the Random Phase Approximation (RPA): 

Xc(Q) = Xo(Q)/[l + V(Q)xo(Q)] with xo(Q) the bare 
susceptibility independent of the Coulomb repulsion and 
V(Q) the interaction potential in reciprocal space at 
wavevector Q. Since the absolute value of V(Q) in the 
model with long-range interactions is much reduced (by 
a factor of 2.5) as compared to the case with nearest- 
neighbor interactions alone, the above argument predicts 
Vco — 4t in the HWM, in agreement with the ED re- 
sult. Remarkably, the three different estimates based on 
ED, PIRG and the RPA scaling argument lead to the 
same value of Vco- This result obtained by consider- 
ing the full long-range Coulomb potential already settles 
a long-standing puzzle in quarter-filled layered organic 
conductors, many of which are found to be metallic even 
for values of V comparable or larger than the bandwidth 
of the material. 

The second critical point corresponds to the metal- 
insulator transition, signaled by the vanishing of the 
Drude weight. The finite-size scaling of this quantity is 
shown in Fig. [5] (b). The Drude weight displays a large 
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drop at about Vmi ~ (7 — S)t after extrapolating to the 
thermodynamic limit. The broad region Vco <V < Vmi 
therefore corresponds to a charge ordered, metallic phase. 
This phase, which can be easily discerned from a metal- 
lic CDW (see below), might be intimately related to the 
hybrid phase predicted at the quantum melting of the 
Wigner crystal in the continuum . 21 ' 22 The robustness of 
the COM phase can indeed be contrasted to the model 
with short range interactions, where a rather narrow in- 
termediate region (if any) arises: between Vco — 
and Vmi — 1.8t also consistent with the COM found with 
cluster dynamical mean-field theory 2 ^. Spin correlations 
are also enhanced in the same parameter range, suggest- 
ing that charge and spin do order together in the metallic 
phase before the metal-insulator transition occurs. Such 
spin order corresponds to antifcrromagnctism between 
charge rich sites due to "ring" exchange processes^. It 
can be noted that using a more realistic value of U = 20t 
significantly increases the value of Vmi, while leaving 
Vco essentially unchanged, which effectively broadens 
the region of stability of the COM phase. 



B. Excitation spectra 

As we proceed to show, the charge ordered metallic 
phase found here shares common aspects of both Mott 
and Wigner physics. Its anomalous properties are a di- 
rect manifestation of the collective nature of its low-lying 
excitations, originating from the long-range Coulomb re- 
pulsion as well as the strong correlation effects associ- 
ated with the on-site Coulomb interaction in a narrow- 
band system. To characterize the excitation spectrum 
we now introduce two representative quantities: the 
single-particle charge gap A c h measured in photoemis- 
sion experiments, that represents the energy required to 
add/remove an electron to/from the system; the optical 
gap A opt defined as the the lowest charge excitation al- 
lowed by the dipolar matrix element, corresponding to 
the onset of finite-frequency absorption in an optical ex- 
periment. The former carries information on the one- 
particle excitations, the latter gives access to the local 
charge fluctuations, pertaining to the collective sector of 
the spectrum. While such quantities coincide for non- 
interacting electrons, their comparison in an interacting 
system provides valuable information on the excitation 
spectrum. These quantities are plotted in Fig. [3Ja) and 
(b) respectively for the EHM and HWM in an L = 16 
cluster, and allow us to analyze the relative role played 
by the single particle against the collective excitations. 

Let us consider the EHM first. At large V t the 
charges order in the checkerboard pattern shown in Fig. 
131 The process of adding an electron or a hole (Fig. 0}t) 
defines the one-particle charge gap A c h = 4V. Reducing 
the interaction strength leads to a progressive closing of 
A c h (the boundary of the yellow region in Fig. [3^,) , that 
eventually triggers the metallization transition and the 
consequent melting of the charge ordered state. In addi- 
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FIG. 3: (Color online) Excitation spectra for the models with 
(a) nearest neighbor (EHM) and (b) long-range interactions 
(HWM) on L — 16 sites. The symbols represent respectively 
the holon (triangles) and exciton (dots) peaks in the opti- 
cal conductivity, and the lowest excitation of finite energy as 
given by the ED (squares). The light shaded (yellow) region 
is delimited by the single-particle charge gap A c h, while the 
dark (blue) region corresponds to the collective excitations 
of the Wigner lattice, bounded from below by A co «. The 
dashed line is the predicted exciton energy in the thermody- 
namic limit, 0.62V (see the appendix). 



tion to the usual one-particle excitations, the motion of a 
particle in the EHM from the perfect checkerboard pat- 
tern to one of its neighboring unoccupied sites gives rise 
to a local (neutral) charge fluctuation depicted in Fig. [4p. 
The latter costs an energy A opt = 3V, which at large V 
lies below the single-particle gap. As V/t decreases, how- 
ever, the quasiparticlc excitations rapidly take over owing 
to their large kinetic energy gain ~ St. As can be seen 
in Fig. in the region around Vmi there is no longer 
a clear-cut distinction between the one-electron and the 
collective sectors, and the neutral charge fluctuations do 
not play a major role in the transition mechanism. In 
fact, the critical point in the case of short range interac- 
tions can be estimated by comparing A c h with the band- 
width w St: this qualitative argument yields Vmi ^ 2t 




FIG. 4: (Color online) Lowest-energy excitations of the 
checkerboard Wigner lattice: (a) Charged (monopoles) , (b) 
neutral (dipolar) and (c) (quadrupolar) defects, (a) cor- 
responds to the process of adding/removing a particle to 
form an interstitial/ vacancy, which determines the charge 
gap A c h — 1.62V; (b) is a local (dipolar) charge fluctua- 
tion of energy A op t = 0.62^, that can be viewed as a tightly 
bound interstitial-vacancy pair, or an exciton; (c) is the next 
(quadrupolar) excitation of energy A qua d = 0.65V . 



(the < sign accounts for the bandwidth reduction due 
to the local U), in agreement with the numerical result 
V M i =i 1-8*. 

This picture changes drastically when the full long- 
range Coulomb repulsion is included, as in the Hubbard- 
Wigner model Eq. (1). The key point is that in this 
case the defects of the checkerboard pattern are no longer 
confined^- Since the electrostatic energy is essentially in- 
sensitive to local details, being mostly determined by the 
long-range tails of the Coulomb potential, several elec- 
tronic configurations can be found which are almost de- 
generate with the local defect of Fig. UJd. For example, 
the energy of the quadrupolar excitation shown in Fig. 
3J; is only 0.03V larger than that of Fig. 5b, making 
these two states practically degenerate for any realistic 
value of V. As a result, the electron can resonate be- 
tween these states, forming a "charge droplet" that can 
propagate coherently at long distances with a net kinetic 
energy gain oc t. This phenomenon is directly reflected 
in the excitation spectrum shown in Fig. [3Jd: a collec- 
tive, itinerant defect state emerges well below both the 
optical gap A opt and the charge gap A c h, whose energy 
A co u ~ A opt — 4.2i qualitatively agrees with the predic- 

A ov t — 2\[2t derived in 



tion of the defect model, A co u 
the appendix. 

We see that due to the derealization of defects in the 
presence of long-range interactions, these collective ex- 
citations can gain a kinetic energy comparable to the 
one-particle excitations. As a result, the separation of 
energy scales between the collective sector and the single- 
particle sector, characteristic of the strongly interacting 
limit, survives down to Vmi and below. Accordingly, the 
metallization transition as well as the resulting charge 
ordered metallic phase in the region Vco < V < Vmi 
appear to be entirely driven by the low-energy collective 
sector—. The finite-size scaling shown in Fig. [5] strongly 
suggests that the one-particle charge gap A c h remains fi- 
nite within the metallic phase: the charge gap is found to 
open concomitantly with the charge ordering transition 
occurring at Vco ~ 4i < Vmi , i- e. in the range of pa- 
rameters where the system has a non-zero Drude weight. 
The region Vco < V < Vmi therefore corresponds to a 
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FIG. 5: Finite-size scaling of the single-particle charge gap, 
A c h, for different cluster sizes. Linear extrapolation to the 
thermodynamic limit indicates that the charge gap opens at 
V ~ At « Vco < Vmi suggesting a phase with a finite charge 
gap, charge order and non-zero Drude weight (see also Fig. 

i). 



charge ordered, metallic phase, with a vanishing single- 
particle weight at low energy. 



C. Spectral probes 

The excitation spectrum described in the preceding 
paragraphs provides a way to univocally characterize the 
anomalous COM phase predicted here. First of all, in 
conjunction with the charge ordering observed by X-ray 
diffraction or NMR, photoemission or tunneling experi- 
ments should find a finite single-particle gap (see Figs. 
[3Jd and[S]), despite a markedly metallic behavior. In fact 
both the large value of the charge gap as well as the 
energy dispersion of the lowest-lying single-particle exci- 
tations are well captured by the defect model presented 
in the appendix. There it is shown that the dispersion of 
the charged interstitials of Fig. is accurately described 
by the following formula: 



E k = E - 2t 1 [cos{2k x ) + cos(2k y )] 

-2t 2 [cos(k x + k y ) + cos(k x - k y )\. 



(2) 



representing the motion of interstitials respectively along 
the principal axes and along the diagonals of the square 
lattice. The physical picture that emerges from the dis- 
persion relation Eq. is that of interstitial defects mov- 
ing as a separate fluid on top of the charge ordered pat- 
tern, which bears a strong resemblance with the hybrid 
phase of the Wigner crystal in the continuum^ An im- 
portant implication of the defect model is that the hop- 
ping parameters t\,t2 governing the motion of intersti- 
tials are essentially constant throughout the charge or- 
dered phase: the one-electron band dispersion depends 
on the interaction strength V only via the rigid energy 
shift Eq, related to the charge gap A c h- This can once 
again be contrasted with the standard CDW picture 
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FIG. 6: The spectral function A(k, w) as obtained from Lanc- 
zos diagonalisation of the Hubbard- Wigner model on a L = 16 
cluster at V/t = 10. The different points in the Brilloun zone 
are T = (0,0), X = (tt.O), A' = (7T,ir). The red dashed line 
is the band dispersion Eq. Q of defects moving along the 
diagonals of the square lattice. The offset energy Eo has been 
adjusted to fit the overall position of the bands. Note that 
the critical values Vco and Vmi at this system size are larger 
than those extrapolated to the thermodynamic limit, so that 
the value V/t = 10 is representative of the charge ordered 
metallic region Vco <V< Vmi. 
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where the band dispersion shrinks as oc t 2 /V upon in- 
creasing V. 

Fig[5] shows the spectral function A(k,oj) obtained 
from the exact diagonalisation of a L = 16 cluster for 

V = lOt, which lies well within the charge ordered metal- 
lic phase for this system size. The main peak clearly 
follows the dispersion of the defect band [Eq. ([2]), with 
t\ = 0, t,2 = — 0.3t] representing the motion of intersti- 
tials along the diagonals of the square lattice. A similar 
agreement applies throughout the charge ordered metal- 
lic phase Vco < V < Vmi , improving at larger values of 
V. 

The calculated optical conductivity also reveals clear 
signatures of the predicted COM phase. As shown in Fig. 
[Tfa) , in addition to a Drudc peak (presumably anoma- 
lous, being possibly carried by collective excitations), the 
optical spectra exhibit two distinct finite-frequency ab- 
sorption bands, confirming the 'mixed' nature of such in- 
termediate phase: a low-frequency 'holon' peak at about 
(3— 4)i originating from the short-range correlations, and 
enhanced by the non-local Coulomb repulsio n 27 i 28 , plus 
a broader excitonic band at A opt = 0.62V, correspond- 
ing to the dipolar defects of the Wigner lattice. Fig. [7jb) 
shows the integrated optical spectra, I(ui), across the CO 
transition. For V = 0, most of the spectral weight is in 
the Drude weight and I{u>) saturates rapidly with oj. As 

V increases a strong redistribution of spectral weight oc- 
curs with the low energy optical weight being transfered 
to higher energies (of the order of A opt ). The total inte- 
grated spectrum is a measure of the kinetic energy of the 



FIG. 7: (Color online) Optical properties across the CO tran- 
sition in the Hubbard-Wigner model, (a) The 'holon' ab- 
sorption band of the strongly correlated metal is strongly 
enhanced by V, dominating the optical conductivity spec- 
tra for V < Vco- Within the charge ordered metallic phase, 
Vco < V < Vmi a broader band at about 0.62V arises due to 
bound excitons (Fig. [4ja) . (b) Integrated spectral weight of 
the optical conductivity for the same parameters as in (a). A 
strong transfer of spectral weight from low energies to high en- 
ergies across the CO transition is evident. Results are shown 
for a L = 16 cluster. 



system as derived from the /-sum rule^. Fig. [TJb) shows 
how the electron kinetic energy is suppressed by V due 
to electron correlation effects close to the CO transition. 



IV. CONCLUDING REMARKS 

The quarter-filled layered organic crystals of the (9-type 
are ideal realizations of a Wigner lattice in which some of 
the existing experiments could be interpreted in light of 
the present scenario. The horizontal charge ordered pat- 
tern observed through X-rays^ below the CO tempera- 
ture T CO = 190K in 0-(ET) 2 RbZn(SCN) 4 , is consistent 
with CO driven by the off-site Coulomb repulsion but 
inconsistent with any Fermi surface nesting vector. In 
the metallic salt 6>-(ET) 2 CsZn(SCN) 4 , both stripe-type 
and three-fold CO patterns are observed, both related 
to strong off-site Coulomb repulsion^. Metallic phases 
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above Tqo display 'bad' metallic behavior with very weak 
T-dependence and absolute values much larger than the 
Mott-Ioffe-Rcgel limil^i indicating mean free paths much 
smaller than the lattice parameter. On application of 
an electric field to 6>-(ET)2CsCo(SCN)4, charge order is 
melted and the conductivity is increased by several orders 
of magnitude 3 ^. 

All the above observations suggest that an unconven- 
tional metallic phase might be realized in such materi- 
als at the melting of the Wigner lattice driven by quan- 
tum fluctuations, which should be further experimentally 
probed. Such a metal should display a charge gap in one- 
electron probes such as photocmission or STM. or a large 
pscudogap on the scale of the non-local Coulomb energy 
V— , together with metallic conduction due to the itiner- 
ancy of defects in the charge ordered background. The 
enhancement of charge fluctuations in the intermediate 
metallic phase could mediate superconductivity 3 ^ close 
to the charge ordering instability occurring, for example, 
in 0-type layered organic conductors and NhSe^. Elu- 
cidating the precise nature of the current-carrying exci- 
tations emerging at zero energy in such anomalous COM 
phase remains an open challenge, as is the understanding 
of the interplay between the Mott and Wigner mecha- 
nisms for the metal-insulator transition^ 
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APPENDIX A: DEFECT DECONFINEMENT IN 
THE WIGNER LATTICE 

A system of electrons is expected to crystallize when 
the mutual Coulomb interactions dominate over the ki- 
netic energy. Minimizing the electrostatic interaction en- 
ergy in a two-dimensional electron gas results in a trian- 
gular arrangement of the charges. The preferred trian- 
gular ordering is altered by the presence of an under- 
lying periodic potential in a solid i 36 ' 37 i 38 for a density 
of n = 1/2 electrons per site on a square host lattice 
the configuration of minimal electrostatic energy is the 
checkerboard order shown in Figf5] We designate such 
"Wigner crystal on an underlying host lattice" simply as 
a Wigner lattice (WL). 

As we show below, the physics of the Wigner lattice 
in the strongly interacting regime V 3> t can be effec- 
tively understood by considering only a small number of 
defects of the classical checkerboard, which play a dom- 
inant role in the low-energy excitation spectrum. These 
defects can be separated in two different classes. The first 
class corresponds to fluctuations of the charge density in- 



duced by quantum fluctuations in the presence of a non- 
vanishing transfer integral t ^ 0, that will be denoted as 
neutral collective excitations. Excitations of the second 
class, termed charged one-particle excitations arise when 
an extra charge is added/removed to/from the system, 
which is the typical situation in a photocmission exper- 
iment. In the following paragraphs we analyze the two 
categories separately. 



1. Neutral, collective excitations 

In the limit of strong Coulomb interactions V ^> t, the 
excitations of the Wigner lattice can be classified in terms 
of the potential (Madelung) energy of the electronic con- 
figurations. The lowest-energy excitations of the checker- 
board pattern are the dipolar and quadrupolar defects 
illustrated in Figf?)} and c. These are obtained by mov- 
ing respectively one or two electrons away from their 
preferred equilibrium position, as shown by the arrows. 
The corresponding electrostatic energies can be evalu- 
ated to arbitrary accuracy through standard Ewald sum- 
mations, being A opt = 0.6155V and A qua d = 0.6453V. 
The next electronic configurations (not shown) have en- 
ergies > 0.898V and will be discarded in the following 
discussion. 

Using the procedure of Ref^, we calculate the en- 
ergy gain induced by a finite transfer integral t between 
molecular sites. This is achieved by restoring the trans- 
lational invariancc and solving the tight-binding problem 
in Fourier space within the reduced subspace consisting 
of the two defects described above. Accounting for the 
two possible orientations of the quadrupolar defect and 
the four orientations of the dipolar defect, we are led to 
diagonalise the following 6x6 matrix (the intermolecular 
spacing a is taken as the unit length): 
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The key point here is that due to the long-range na- 
ture of the electron-electron interactions, the dipolar 
and quadrupolar defects are almost degenerate in en- 
ergy: A qua d — A opt — 0.03V. As a result a resonant 
tunneling is established between the two non-equivalent 
states for any physically reasonable value of the V/t ra- 
tio. The electron density is then shared between these 
states, forming an "electron droplet" . Such delocalisa- 
tion of the charge yields a kinetic energy gain cx t and 
allows the droplet to propagate coherently at long dis- 
tances with a band-like dispersion. The present scenario 
can be contrasted to the model with nearest-neighbour 
interactions, where the defects remain strongly localized 
in space unless the transfer integral overcomes the large 
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FIG. 8: (Color online) Comparison of the optical conductivity 
calculated at V = 17i from ED of a L = 16 cluster with the 
one obtained from the defect model, including respectively 
the 6 and 12 lowest-lying states. 



energy barrier separating the states depicted in Figs. [4)d 
and c, i.e. t > {A quad - A opt ) = V. 

By setting A qua d = A opt in the matrix (|A1[) we obtain 
an analytical expression for the energy gain of the droplet 
state at k = 0: 

A fc=0 = A opt - 2V2t. (A2) 

This form is in satisfactory agreement with the exact di- 
agonalization (ED) result: Ak=a — A opt — 4.2*, indicat- 
ing that the physical mechanism of defect derealization 
is correctly captured already in the small droplet approx- 
imation, where only dipolar and quadrupolar states are 
included. 

We note that while the quadrupolar states are the ones 
responsible for the band dispersion, it can be shown fol- 
lowing the lines of Refi^ that the optical spectral weight 
is mostly carried by the dispcrsionlcss dipolar excita- 
tions. As a result, the peak in the optical conductivity 
remains centered at an energy uj = A opt while the peak 
at to = Afc = o is not optically active. The optical con- 
ductivity calculated from the defect model at V = 17* 
is shown in Fig. [5J and compared with the spectrum 
obtained from the full ED in a L = 16 cluster. The mini- 
mal model of Eq. (|A1|) including 6 defect states is able to 
reproduce the correct position of the exciton peak. Ex- 
tending the defect subspace to a total of 12 states also 
accounts for the emergence of a sideband below the main 
peak, as seen in the numerical data. 



2. Charged, single particle excitations: interstitials 
and vacancies 

An analogous procedure can be carried out to deter- 
mine the dispersion of the one-electron states that are 
probed in photocmission experiments. It is convenient to 
include a static compensating charge density of +e/2 per 




FIG. 9: The lowest-lying states with one electron added to 
the perfect checkerboard pattern. From left to right the cor- 
responding Madelung energies in an infinite system are (a) 
Ei = 0.80781/, (b) E D = 0.9233 and (c) E M = 1.1304. 



site in the calculation (equivalent to the usual compen- 
sating jellium in continuum models), which restores the 
overall charge neutrality and the particle-hole symmetry. 
We can therefore focus on the addition spectrum alone, 
since the removal spectrum is obtained by symmetry. 

As usual we start from the V 3> * limit. Adding an 
electron to an empty site of the Wigner lattice creates an 
interstitial, of energy Ej = 0.8078V. The next low-lying 
states with N + 1 electrons have energies respectively 
E D = 0.9233 and E M = 1.1304. All these states are 
illustrated in FigJSJ The charge gap is defined as A c h = 
E(N+1) + E(N-1)-2E(N) which gives A ch = 2\Ei\ = 
1.61551/. 

In order to study the effect of quantum fluctuations 
induced by a finite transfer integral * we evaluate the ac- 
tion of the kinetic energy operator within the subspace 
of the defect states defined above. Considering the mul- 
tiplicities arising from the possible orientations of these 
states (respectively 2,4,8) one obtains the 14 x 14 ma- 
trix shown in Table IB 21 The charge gap at finite * is 
obtained from the energy of the lowest state at k = 0. 
In the region 10 < V/t < 20 it follows the linear form 
A c h — 1.8V — 9.6*, in good agreement with the ED value: 
A ch ~ 1.8V -7.0*. 

To get more insight on the one-particle excitations of 
the Wigner lattice, we observe that the bands arising 
from the diagonalization of the matrix in Table IB 21 are 
accurately described by the following formula: 

E k = E Q - 2fi[cos(2k E ) +cos(2k y )] (A3) 
— 2*2[cos(fc.j; + ky) + cos(k x — k y )]. 

An important implication of the defect model is that the 
hopping parameters *i,*2 governing the motion of de- 
fects are essentially constant throughout the charge or- 
dered phase: the one-electron band dispersion depends 
on the interaction strength V only via the rigid energy 
shift Eq. For a direct validation of the defect model, 
FigJB] shows how the spectral function A(k, oj) obtained 
from the exact diagonalisation of a L = 16 cluster within 
the charge ordered metallic phase closely follows the dis- 
persion of the main defect band [Eq. (|A3|) . with *i = 0, 
*2 = —0.3*], which represents the motion of interstitials 
along the diagonals of the square lattice. 
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APPENDIX B: FINITE-SIZE EFFECTS 

As was pointed out in Ref.— , in addition to the usual 
size effects encountered in numerical studies of model 
systems on finite clusters (hereafter denoted as quantum 
finite-size effects) , specific issues arise when dealing with 
the long-range Coulomb potential. These originate from 
the inaccuracy of calculating the electrostatic energy of 
the electronic configurations in a finite system. Such ef- 
fects are termed classical finite size effects and become 
relevant within the charge ordered phase at large V/t. 

1. Quantum finite-size effects 

Ordinary finite size effects in the EHM have been 
explored previously by comparing results of different 
quantities with different number of cluster sites L = 
8, 10, 16 and 20^. Other possible ways of evaluating fi- 
nite size effects are based on averaging over the bound- 
ary condition^. In the case of the extended Hubbard 
model both procedures are found to lead to similar crit- 
ical values^. 

We have performed finite-size scaling for the HWM 
by evaluating the charge correlation function, C(7r,7r), 
Drude weight, D, and charge gap, A c /, on L = 
8, 10, 16, 18 clusters, and analyzing the extrapolation to 
the thermodynamic limit. In Fig. [5] we show C(7r, 7r) 
plotted as a function of 1/L and the Drude weight, D, 
dependence with l/VT. The finite size scaling of the 
single particle charge gap, A c ^, is shown in Fig. [5] Our 
finite-size scaling analysis strongly suggests the presence 
of a broad metallic phase with both charge order and a 
finite single-particle charge gap for Vco <V< Vmi- 



amounts to performing an infinite periodic repetition of 
the finite simulation cell. Although a system of infinite 
size is effectively recovered through this procedure, this 
gives rise to spurious contributions to the energy arising 
from the interaction of a given electronic configuration 
with its images on the repeated cells. Since the simula- 
tion cell is overall neutral, the interaction between equiv- 
alent cells has a dipolar nature and the corresponding 
error in the energy scales at most as AE ~ V/L 3 ^ 2 (L 
being the total number of sites). 

Having identified in Appendix [X] those electronic 
configurations that are mainly responsible for the low- 
energy behaviour of the Wigner lattice, we can precisely 
address the effect of such uncertainties on the physics 
of interest here. To this aim we report in Table U the 
electrostatic energies of the configurations used in the 
defect model of Appendix [XJ and compare them with 
the corresponding values in a 4 x 4 cluster (energies are 
expressed in units of V) . As we can see, finite-size errors 
on these states are typically of the order of 5 — 10% 
or less, which gives us good confidence on the final 
numerical results. 





E op t E qU ad Ej Ed Em 


extended 
4x4 


0.6155 0.6453 0.8078 0.9233 1.1304 
0.5776 0.6258 0.8078 1.0370 1.1207 



TABLE I: Electrostatic energies of the relevant defect states. 



2. Classical finite-size effects 



Classical finite-size effects are minimized by evaluat- 
ing electrostatic energies with Ewald summations, which 
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TABLE II: The 14 x 14 matrix used to calculate the one-particle excitation spectrum in the defect subspace. 



Ei 





— 2t cos kx 


— It cos ky 










. &x ~\~ 


. k x ~\~kv 


k +k 


k -k 


k -k 


k -k 


k -k 


—te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 





Ei 


— It cos k x 


—It COS ky 
























-te 1 2 


-te 1 2 


-te % 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


-te 1 2 


—It cos k x 





E D 



































— 2t COS fty 


U 


A 

u 


&D 


n 
u 


A 

u 


n 
u 


n 
U 


A 

u 


A 

u 


n 
u 


n 
U 


n 
U 


n 
U 


o 


— 2£ cos k x 








Ed 
































—2t COS fey 











Ed 


























. k x + k v 


■ &x ~\~ 


























-te 1 2 


-te 1 2 














Em 























■ k-x 


■ "-a; ~\~ 


























-te 1 2 


-te 1 2 

















E M 




















■ "-an "I - ™!/ 


■ &x ~\~ 


























-te 1 ^^ 


-te 1 2 




















Em 



















■ i^x m ^~> a y 


























-te 1 2 


-te 1 2 























Em 














. f^x — 




























-te 1 2 


-te 1 2 y 


























Em 







































-te 1 — 5-^- 


-te 1 ^T JL 





























Em 




































-te 1 2 " 


-te l 2 V 
































Em 





. &x — "■y 




























-te 1 2 


-te 1 2 v 



































Em 



\ 



